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Abstract 



_^ We have investigated the fluorescence emission in the Near Infrared from the air and its main components, nitrogen and oxygen. 
^-^j The gas was excited by a 95 kV electron beam and the fluorescence light detected by an InGaAs photodiode, sensitive down to 
about 1700 nm. We have recorded the emission spectra by means of a Fourier Transform Infrared spectrometer. The light yield was 
£N| also measured by comparing the Near Infrared signal with the known Ultraviolet fluorescence, detected by a Si photodiode. 

The possibility of using the Near Infrared fluorescence of the atmosphere to detect Ultra-High-Energy Cosmic Rays is discussed, 
(1) showing the pros and the cons of this novel method. 
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1. Introduction 

The fluorescence of the nitrogen gas in the atmosphere in the 
visible (VIS) and ultraviolet (UV) regions was independently 
proposed in the sixties of last century by Greisen [1 1, Delvaille 
et al. [2|, Suga [3|, and Chudakov [4| as a way to detect high 
energy showers created by cosmic rays impinging on the earth 
atmosphere. Nowadays the air fluorescence method is a well 
established technique used by the experiments HiRes [5 1, which 
is the continuation of the pioneering Fly's Eye experiment |6|, 
AUGER |7|, and by the Telescope Array project |8|. 

The main emission of the fluorescence light is in the UV re- 
gion between 300 and 450 nm. The transparency of the atmo- 
sphere at those wavelength is limited by the presence of molec- 
ular oxygen and ozone, which absorb the UV photons. The 
Rayleigh scattering, whose cross section scales as 1/A 4 , and the 
Mie scattering further reduce the number of photons which ar- 
rive at the detector. We introduce an extinction length A(A), 
function of the wavelength A, which enters in the transmission 
of the light through the relationship I(x, A) = Iq ■ exp(-x/A{A)), 
where I(x, A) is the intensity of a collimated light beam of wave- 
length A after having travelled for the distance x in the atmo- 
sphere, and Iq is the initial intensity (at x = 0). For the UV, 
A < 10 km. 

This fact poses two main drawbacks. First, the detector can- 
not be placed at sea level, but at altitude. For example, AUGER 
is at about 1400 m above sea level, HiRes at 1700 m. 

Second, the detection range of the instruments extends to a 
distance of the order of A. AUGER, for example, is capable to 
detect the highest energy air showers at a distance of the order 
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of 30 km |9|. This is correlated to the total observation rate, 
which goes roughly as A 2 , which is therefore limited also by 
the atmospheric transmission. 

The atmosphere is very transparent in some regions of the 
near-infrared (NIR) region, below 3 ymi. The molecules which 
present absorption bands are water and, in a minor extent, CO2. 
Several high transparency windows are present: from 0.8 to 
0.9 pmi, from 0.95 to 1.05 ymx, from 1.15 to 1.3 fim, from 1.5 to 
1.8 fim, from 2.0 to 2.4 fxm (see, for example, ref. lfTOl ). More- 
over, the Rayleigh scattering is negligible in the NIR, because 
of the I /A 4 dependence. If the air fluoresces at those wave- 
lengths, then there is the possibility to detect high energy cos- 
mic rays impinging on the atmosphere at a very long distance 
from the detector, thus increasing the observation rate. 

The spectroscopy of the molecule of nitrogen is rather com- 
plex and has been extensively studied in the last century. An 
exhaustive review can be found in ref Jlli Despite of this, to 
our knowledge there are no investigations on the infrared flu- 
orescence emission from N2 or air excited by ionizing radia- 
tion. The only research which extends somewhere in the NIR 
region was performed in the sixties by Davidson and O'Neil 
[12 1, who used photo multipliers to reach the maximum wave- 
length of 1050 nm. 

In this paper we present our first measurements of the fluores- 
cence light in the NIR region of air, nitrogen, and oxygen. The 
light emitted by the gas, excited by a 95 kV electron beam, was 
detected by an InGaAs photodiodes, sensitive down to about 
1.7 fj.m. Spectra have been recorded at atmospheric pressure 
and room temperature by mean of a Fourier Transform Infrared 
Spectrometer. We also report on the light yield, which has been 
obtained by comparing the NIR signal with the known fluores- 
cence output in the UV region, measured by means of a silicon 
photodiode. 
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Figure 1: A schematic view of the experimental setup (not to scale). Not shown are the vacuum and the gas lines. 



2. Experimental details 

A schematic of the experimental setup is shown in figfl] 
An electrostatic gun 1 13] accelerated electrons, emitted by a 
hot filament, up to 100 keV kinetic energy. Electrons exited the 
gun through a CVD diamond window [ 14], 20 pm thick, whose 
external side had been metalized to screen the fluorescence light 
emitted by the diamond itself. The gun was operated at 95 kV 
in pulsed mode, with pulses 300 ps long and rate of 50 Hz. 
In those conditions, the electron current behind the diamond 
windows was a few hundreds pA. 

The chamber containing the gas was a stainless-steel cylinder 
with the axis along the beam path and volume about 300 cm? . 
The chamber size was a compromise between the range of the 
electrons in the gas and the need of an adequate solid angle for 
the detection of the light. The ionization cloud was almost en- 
tirely contained in the chamber, and the electrons that hit the 
materials of the chambers did not affected neither disturbed the 
measurements, as proved in section 3. Perpendicular to the 
axis, one opposed to the other, two quartz windows allowed 
to view inside the chamber. They were 10 mm thick in order to 
stop the X-rays generated in the chamber, that could alter the 
signals of the photodiodes. By construction, they looked ex- 
actly at the same gas region, a small fraction (about 5%) of the 
chamber volume. The light to the spectrometer exited through 
a third quartz window on the closing flange of the vessel. In- 
side the chamber, an electrode collected the ionization charge 
for the beam monitoring and normalization, which was fed into 
a charge amplifier. The output signal Qbeam was displayed on a 
digital oscilloscope and also used as a feedback to stabilize the 
electron gun current against long term drifts. Qbeam was com- 
pared to a reference voltage level, and the difference signal was 
used to opportunely bias the grid in front of the filament, which 
regulated the electron beam current intensity before accelera- 
tion. 

Gas from certified bottles was flown in the chamber at room 
temperature and atmospheric pressure with a flux of 0.2 i/min. 
We used N2 (impurity concentration < 5.5 ppm), O2 (impurity 
concentration < 5 ppm), and dry air, which is a mixture of 80% 



nitrogen and 20% oxygen (impurity concentration < 10 ppm). 

The fluorescence light was detected by two solid state pho- 
todiodes (PDs), which could be placed on the two lateral win- 
dows of the chamber or inside the spectrometer. The UV light 
was detected by a squared Si PD 1 15 1, 10 x 10 mm 2 area, with 
a quantum efficiency (QE) of about 55% from 300 to 400 nm. 
The infrared radiations was detected by a round InGaAs PD 
[ 16 1, 5 mm diameter, sensitive from about 500 nm to 1.7 pm 
and cooled down to about -40 °C by a built-in Peltier cooler. 
Its QE is above 80% from 1.1 to 1.5 pm. The QE curves of 
both PDs are shown in figj2] Each photodiode was readout by 
a dedicated charge preamplifier, based on the UA1 hybrid ifTTl . 
The gain was measured by injecting a known calibration charge. 
The output signals were displayed and measured by means of a 
digital oscilloscope. 

Since the Si PD is also sensitive to the NIR radiation, we 
used some colored glass filters to select only the UV photons. 
The combinations of two filters, the filter FGB25 (by Thorlabs) 
and the filter HA30 (by Hoy a Corp.), produced the transmission 
curve shown in fig|2] Only two bands are allowed, from 300 to 
500 nm, which is the region of the UV fluorescence, and from 
700 to 900 nm, where the fluorescence light output of air and 
nitrogen is negligible |[T2l . 

Spectra were recorded and analyzed by a Fourier Transform 
Infrared (FTIR) spectrometer (mod. Equinox55 by Bruker Op- 
tics) which is essentially a two-arm Michelson interferometer 
with a 632 nm He-Ne laser as reference light source. For the 
Equinox55 the minimum measurable wavelength was 632 nm. 
For details of this well known spectroscopic technique see for 
example ref.[18 19]. The spectrometer was controlled and the 
data were elaborated by the software OPUS 5.5 (by Bruker Op- 
tics). 

The laser light covered the same path of the radiation under 
measurement, and reached the detector, producing a very large 
noise. To filter the He-Ne light, we placed a colored glass filter 
in front of the photodiode, which transmitted only photons with 
wavelength longer than 780 nm (filter FGL780 by Thorlabs). 

The FTIR spectrometer recorded the interferogram of the 
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Figure 2: Quantum efficiency curves of the PDs. Dashed line = Si PD; full line = InGaAs PD. It is also shown the total transmission curve of the colored filters used 
to isolate the UV region. 



light source as a function of the wavenumber k - l/A. In 
the calculation of the Fourier transform of the interferogram 
to obtain the spectrum, we applied the apodization function 
Blackman-Harris (3-terms), and the phase correction with the 
Mertz method (no peak search) lfT8l[T9l . 

The spectrometer resolution on k ranged between 20 and 
5 cmT x but the final spectrum resolution was worse because of 
the low spatial coherence of the light source and of the math- 
ematical operations involved in the interferogram-to-spectrum 
conversion. Note that a constant resolution 6k on k implies a 
non constant wavelength resolution SA, given by 6 A = A 2 ■ 6k. 

To reduce fluctuations and noise in the spectrum, the final re- 
sult was obtained by averaging together many different spectra 
(between 10 and 20). 

3. Results 

As preliminary operation, we verified that no light emission 
occurred from the CVD diamond window or from the quartz 
windows or other materials of the chamber. To accomplish that, 
we evacuated the chamber with a rotative pump and run the 
electron gun at the highest beam current and energy. We could 
not observe any light signal. 



3.1. Spectra 

Figures [3] [5] and [6] show the measured spectra of nitrogen, 
oxygen and dry air, respectively. The spectra were recorded 
at atmospheric pressure and room temperature (298 K), and 
measured with a resolution 6k — 5 cmr 1 , apart from oxy- 
gen. They are corrected for the transmission curve of the filter 
FGL780 and for the QE curve of the InGaAs PD, and truncated 
to 1650 nm. To identify the atomic transitions, we used the 
NIST Atomic Spectra Database Il20l . 

3.1.1. Nitrogen 

Molecular nitrogen has a very rich spectrum, with numerous 
band systems corresponding to electronic transitions of N2 and 
Nj- The near infrared region is dominated by the First Posi- 
tive System (B 3 IL g - A 3 Z^) (FPS). A less important system is 
the A 2 n„ - X 2 I,g Meinel System (MS) of the ionized nitrogen. 
The Herman Infrared System (HIS) also appears, which extends 
from 700 to 910 nm ifTTl . 

The spectrum in figj3] is dominated by the unresolved band 
around 1040 4- 1050 nm, which has been already observed by 
Ifl2l . The band is the overlap of two nearby atomic transitions 
at 1040 nm, characteristic of the forbidden doublet transition 
[N /] 32 , and of the FPS (0-0) band (figg). The second more 



intense band is the (0-1) band of the FPS centered at 1230 nm. 
Other minor bumps has been identified as belonging to the HIR, 
to the MS (at about 1187 nm) and to the FPS (the (0-2) band 
centered at 1485 nm). Two very weak atomic lines of the neu- 
tral nitrogen atom are distinguishable at the side of the main 
band, at 101 1 nm (transition 3d 4 F — > 3p 4 D°) and at 1076 nm 
(transition 3d 4 P -> 3p 4 P°). 

Because of the weak intensity and the poor resolution, bands 
from about 1250 to 1450 nm cannot be unambiguously identi- 
fied. 

3.1.2. Oxygen 

Since the oxygen fluorescence light output is sensibly lower 
than the nitrogen or dry air one, we run the spectrometer with a 
resolution 5k = 20 cm -1 . 

Molecular oxygen is a weak light emitter since transitions 
from excited states to the ground state are strongly forbidden, 
for most of the excited states. As a fact, the only line we could 
clearly identify to belong to a molecular transition is the faint 
line at 1268 nm, which is the well-known A g — > 3 E g magnetic 
dipole intercombination transition ETTl . All other lines have 
been identified as atomic transitions from the excited neutral 
oxygen (O I), created by electron-impact dissociation of the O2 
molecule, or from the single ionized atom (O II). 

The list of the lines and relative transitions is reported in ta- 
ble [T] The line at 1130 nm is the sum of two very close lines, 
which were not resolved by the spectrometer. The line inten- 
sities, after the continuum subtraction, are normalized to the 
highest intensity line. 



and V[r, respectively. Preliminarily, the symmetry of the appa- 
ratus was verified by exchanging the two detectors and check- 
ing that their signal amplitude did not change. The data were 
taken at different beam intensities to verify whether saturation 
effects or other effects due to the high ionization density of the 
gas, are present. The beam intensity was measured through the 
charge Qbeam collected by the beam pickup, which is propor- 
tional to the energy deposited in the gas. Each point was the 
average of about a hundred signals, performed by the digital 
oscilloscope. As can be seen from fig[8] the dependence of Vuv 
and Vir on Qbeam is proportional to the charge over 3 orders 
of magnitude of beam intensity, which consequently excludes 
non-linear effects due to the ionization density. The photodi- 
ode signals had a small contribution from the X-rays generated 
inside the chamber and not absorbed by the quartz windows. 
To quantify the effect, we placed in front of the photodiodes a 
thin black paper foil, which is opaque to the NIR light but com- 
pletely transparent to the X-rays, and measured the signal vs. 
Qbeam- The X-ray signal (fig[8]> is proportional to Qbeam m the 
whole explored range. The effect is however small, being about 
2% for SiPD and 17% for InGaAs. We can therefore write 



Va(Qbeam) - A a ■ Qbeam + B a • Qbeam (<* - UV, IR) 

where the B a coefficients are related to the X-ray signals and are 
obtained from the fits of fig|8] The A a coefficients are related 
to the total photon yield Yir and Yuv in the near-infrared and in 
the UV region through the relationships: 



3.1.3. Dry air 

From the comparison with the spectra of N2 and O2 (fig|3]and 
fig|6]and also the expanded views fig|4]and fig|7]), it is evident 
that the spectrum of the dry air is not simply the weighted sum 
of the two. Atomic and molecular interactions among differ- 
ent elements change the population of the excited states of the 
species and the de-excitation spectrum. The spectrum of dry 
air is densely populated with bumps and peaks, which could 
not always be easily identified. It is dominated by the nitro- 
gen FPS band (0-0) at 1046 nm, but, contrary to pure N2, the 
doublet atomic lines at 1040 nm are absent (compare fig|4]and 
fig|7]i. Many atomic lines of the excited neutral nitrogen atom 
are present, some of which were not observed in the pure gas, 
or have a different intensity, such as, for example, the two lines 
at 1 100 nm and 1076 nm at the sides of the FPS(O-O) band (fig|4] 
and [7}. The only lines ascribable to oxygen are the O I lines at 
845 and 1130 nm. Many bumps are probably bands of the N2 
FPS. Since an exhaustive discussion of all the features of the 
spectrum is beyond the purpose of this paper, we do not go into 
further detail. In table [2] we report the list of the main peaks 
together with their associated transition (when identified), and 
their intensity. 

3.2. Light yield 

To measure the light yield Yir of the dry air in the NIR region 
we compared the signal output of the Si and InGaAs PDs, Vuv 
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C • Yuv • &uv ■ GsiPD ' £uv 



Air - C ■ Y IR ■ Q. IR ■ Gj„CaAs ' ^IR 



(1) 



(2) 



where C is a constant, Cluv and Q.ir are the solid angles for the 
Si and InGaAs detectors, respectively, GsiPD and G/„caAs the 
gains of the two charge amplifiers, is the average quantum 
efficiency e a (A) of the photodiode (fig[2j» weighted by the ap- 
propriate light fluorescence spectrum S a (A), which must also 
include the transmission curve T a (A) of the optical filter when 
present: 



/ e a (A)-T a (A)-S a (A)dA 
J S a (A)dA 



(a = UV, IR) 



with Tjr(A) = 1 and T UV (A) given in fig|2] S jr(A) is the ex- 
perimental spectrum of fig|6]while S uvW has been taken from 
(22). 

From the ratio of eq.([T]) and (j2j), we obtain Y IR as: 

„ „ A;r Qj/y GsiPD ^UV 



'IR 



l UV • 



A-UV &IR G InGaAs £ IR 

The quantities Air and Auv have been obtained from the linear 
fits of fig[8] after X-ray background subtraction. The result is 
0.231 ±0.002. Since: 



A 

Auv 



■■ 4.19 ±0.15 



Table 1 : Observed line in the oxygen spectrum. O I is the neutral oxygen atom, O II is the single ionized atom. For the calculation of the intensity, the continuum 
has been subtracted. 



wavelength (nm) transition species intensity 

822 3p 3 D X2 -> 3s' 3 D" OI 10 

845 3p 3 P^>3s 3 S° OI 24 

927 3d 5 D°^3p 5 P OI 6 

1130 3d 3 D° -» 3/? 3 P(1129nm) OI 100 

4s 5 S Q -» 3/? 5 P(1130 nm) 

1246 5/ 5 F ->3d 5 D° O I 5 

1268 1 Ag Z g 2 7 

1301 4p 2 D°^>4s 2 P Oil 13 

1316 4s 3 S°^>3p 3 P OI 29 



Table 2: Observed peaks in the dry air spectrum. When identified, the transition and the atomic/molecular species are reported. For the calculation of the intensity, 
the continuum has been subtracted. 



wavelength (nm) transition species intensity 
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Figure 3: Fluorescence spectrum of nitrogen at room temperature and atmospheric pressure. The bands and lines clearly identified are shown: HIS = N2 Herman 
infrared system; N I = atomic neutral atom; FPS = N2 first positive system; MS = Meinel system. 



= 0.500 ± 0.006 
euv = 0.35 
Wr = 0.80 
the ratio of the light yields is: 

— = 0.21 ±0.03 

Yuv 

The error in the electronic gain ratio is dominated by the uncer- 
tainty of the calibration charges. Since the QE curves of the Si 
and InGaAs PDs are rather fiat in the region of our interest, the 
errors on ejjv and i/R are small and can be neglected. 

The value of Yjjv has been calculated from the measurements 
available in literature, performed in the wavelength region from 
300 to about 430 nm at atmospheric pressure and room temper- 
ature. The data considered by us and their references are listed 
in table [3] Regarding the value in l23l . the authors reported the 
yield as number of photons per unit length. To excite the gas, 
they used a 90 S r /3~ source, which emits electrons with aver- 
age energy of 0.85 MeV. Such an energy deposit is equivalent 
to 0.2059 MeV/m [24] . Then we can calculate the yield per 
deposited energy as reported in the table. 

We calculated the weighted average, which yields Yjjy — 
19.88 ± 0.51 ph/MeV. Therefore, our estimate of the NIR light 



yield Y tR is 

Y IR = 4.17 ±0.53 ph/MeV 

4. Discussion 

The light yield of dry air in the NIR is about 1/5 of the UV 
yield. The QE of the InGaAs semiconductor is about 4 times 
higher than the UV photocathode of a photomultiplier, so the 
total number of photoelectrons is approximately the same in the 
two cases. However, from a practical point of view, when com- 
pared to photomultipliers, operation with InGaAs photodiodes 
poses several drawbacks, such as: 

- small dimension (area < 1 cm 2 ); 

- no multiplication; 

- need of a low noise electronics. 

Therefore, detection of cosmic showers with InGaAs PDs is 
nowadays not compelling. 

The situation could be different if we adopt Si semiconductor, 
whose QE can extend down to the neighborhood of 1100 nm, 
being limited by the l.leV band-gap of the ordinary bulk sili- 
con. Moreover, photodiodes can be produced with an internal 
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Figure 4: Expanded view of the fluorescence spectrum of nitrogen at room temperature and atmospheric pressure in the wavelength region 1000 4- 1100 nm. The 
spectrum is characterized by the overlap of two nearby atomic transitions at 1040 nm (not resolved), and of the FPS (0-0) band. 



gain, obtaining Avalanche Photodiodes (APDs), with gain in 
the range 10 2 -r 10 4 and surface area larger than 1 cm 2 [27|. 
However, limiting the wavelength range would decrease the 
light output. The amount of NIR light with A < 1100 nm is 
about 35% of the total light yield in the 800 -f 1650 nm range. 
Therefore, Si APD in the NIR is about a factor of 10 less sensi- 
tive with respect to the UV light, at parity of detecting surface 
area. 

So far we have not mentioned the main argument in favor 
of the NIR fluorescence, namely, the transparency of the atmo- 
sphere. The UV light produced by a shower 20 km far away 
from the detector is attenuated by the atmosphere by about a 
factor of 10. Then, since atmosphere is transparent in the NIR, 
detection of the NIR fluorescence with a Si APD is, in this case, 
equivalent to the detection of the UV light with photomultipli- 
ers. 

For the general case, it is likely that the detection of UHE- 
CRs through the NIR fluorescence with Si APDs is not much 
disadvantaged with respect to the UV fluorescence with pho- 
tomultipliers. The principal pro for the physics of UHECRs 
would be the increase of the event rate, which scales as A 2 . 

On the opposite side, however, it must be demonstrated that 
the large background noise of the atmosphere in the NIR, which 



is higher than in the UV 11281 . can be handled and does not 
introduce other limitations. 

The technological research on Si APDs towards the enhance- 
ment of the response in the NIR is progressing and has already 
reached remarkable results, not yet transferred on commercial 
devices. For example, several groups have obtained APDs with 
QE above 50% for A > ljjm by laser treatment of the Si surface 
ll29l [30l [3TI . Others used inclusions of SiGe absorbing layers 
to enhance the QE well beyond the Si cutoff wavelength [32|. 
It is therefore thinkable that in the near future the Si APDs are 
much more performing than the actual ones. 

For sake of completeness, we have also to mention the pos- 
sibility to use NIR photomultipliers, whose sensitivity now ex- 
tends down to 1700 nm l33l . They have a small active area 
(smaller than the InGaAs PD) and a QE < 1%. Therefore op- 
eration with such photomultipliers seems not competitive with 
respect to large area Si APDs. 

5. Conclusions 

We have measured for the first time the fluorescence light 
of the air in the near infrared region. The NIR photon yield 
is 1/5 of the UV yield. We propose such fluorescence light as 
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Figure 5: Fluorescence spectrum of oxygen at room temperature and atmospheric pressure. Arrows show the main lines reported in table[T] 



a new way to detect UHECRs. Since the atmosphere is more 
transparent in the NIR than in the UV region, there is reasonable 
expectation that the new method could lead to an increase of 
the UHECR observable rate. This paper is the first step of a 
research program which aims to this goal. 
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Figure 6: Fluorescence spectrum of dry air at room temperature and atmospheric pressure. 




Table 3: Light yields reported in literature and used to calculate Yuv- 

wavelength range (nm) UV light yield (ph/MeV) reference 



300 - 428 
300 - 430 
300 - 420 
290 - 440 



19.67 ± 0.68 
20.38 ± 0.98 
20.8+ 1.6 
17.6 ±2.3 



1 23 1, see text 
[24| 

ma 



9 



FPS(O-l) 




1000 1020 1040 1060 1080 1100 

wavelength (nm) 

Figure 7: Expanded view of the fluorescence spectrum of dry air at room temperature and atmospheric pressure in the wavelength region 1000 ~ 1100 nm. The 
spectrum is characterized by the FPS (0-0) band and by two atomic transitions of the neutral nitrogen. 
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Figure 8: Dependence of the light signal of the Si and InGaAs PDs, Vjjv and V :R respectively, upon the charge collected by the beam pickup. Circles = Si PD data; 
squares = InGaAs PD data; Full symbols = total PD signals; open symbols = X-ray contribution. The full lines are the least square fits to the data with the y = ax 
function. 
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